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ABSTRACT 

We report extensive spectroscopic and differential photometric BVRI observations of the active, 
detached, 1.309-day double-lined eclipsing binary IM Vir, composed of a G7-type primary and a K7 
secondary. With these observations we derive accurate absolute masses and radii of M\ = 0.981 ± 
0.012 Mo, M 2 = 0.6644 ± 0.0048 M 0) i?i = 1.061 ± 0.016 R Q , and R 2 = 0.681 ± 0.013 R Q for the 
primary and secondary, with relative errors under 2%. The effective temperatures are 5570 ± 100 K 
and 4250 ± 130 K. The significant difference in mass makes this a favorable case for comparison with 
stellar evolution theory. We find that both stars are larger than the models predict, by 3.7% for the 
primary and 7.5% for the secondary, as well as cooler than expected, by 100 K and 150 K, respectively. 
These discrepancies are in line with previously reported differences in low-mass stars, and are believed 
to be caused by chromospheric activity, which is not accounted for in current models. The effect is not 
confined to low-mass stars: the rapidly-rotating primary of IM Vir joins the growing list of objects 
of near-solar mass (but still with convective envelopes) that show similar anomalies. The comparison 
with the models suggests an age of 2.4 Gyr for the system, and a metallicity of [Fe/H] «s —0.3 that is 
consistent with other indications, but requires confirmation. 

Subject headings: binaries: eclipsing - binaries: spectroscopic - stars: late-type - stars: fundamental 
parameters - stars: individual: IM Vir 



1. INTRODUCTION 

Our knowledge of stellar structure and evolution rests 
heavily on the comparison between theory and obser- 
vation. Double-lined eclipsing binaries (hereafter EBs) 
have long been at the center of this process, since they 
allow the mass - the most fundamental of all stellar prop- 
erties - as well as the radius to be determined to very 
high precision (and accuracy), often as good as 1%, inde- 
pendently of the distance and independently of any cal- 
ibrations. Such high precisi on enables string ent tests of 
theory, as described, e.g., by I Andersen! (|1991f ). In the last 
decade or so it has become clear that stars in the lower 
main-sequence show significant discrepancies when com- 
pared to standard models. Studies of several key systems 
have shown unambiguously that the radii predicted by 
the models are systematically up to ~ 10% too small, and 
the temperatures ~5% too high (see iRibas et al.l [20081 
and references therein). The deviations are commonly 
attributed to the high level of chromospheric activity 
present in these systems. Orbital periods are typically 
short, and as a result tidal forces tend to synchronize the 
components' rotation with the orbital motion. The rapid 
rotation, in turn, induces the activity, which can mani- 
fest itself in the form of copious X-ray emission, flaring, 
Ha and Ca II H and K emission, spottedness, and other 
effects. 

According to the recent compilation by iTorres et al.l 
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(2009), among the eclipsing binaries with at least one 
component below 0.8 M only five have reliable mass 
and radius determinations with relative errors that are 
smaller than 3%. These are, in o rder of decreasing mass, 
UV P sc (|Po ppeiill997h. GU Boo (lLopez-Morales fc Ribasl 
2005J), YY Gem ([Torres fe Ribasl 12002ft. C U Cnc (IRibas! 
2003), and CM Dra ([Morales et all 12001) . A number 
of other similar systems are known, but are not yet at 
the same level of precision and accuracy. For further 
progress it is therefore important to either improve the 
determinations in the latter cases, or to find new ones 
where those conditions are met. The eclipsing system 
reported here is one of these, which in addition presents 
the largest difference in mass between the components, 
providing in principle extra leverage for the comparison 
with models. 

IM Virginis (also HD 111487, IE 1247.0-0548, a = 
12 h 49 m 38. s 70, S = -6°04'44."9, J2000.0; G7 V, V = 
9.57) was detec ted as an X-ray source with the Einstein 
Observatory by iHelfand fe CaillaultJ (IT982T). Th e radial- 
velocity variability was found by iSilva et aTl (|1987h . 
and subsequen t spectroscopic and photome tric studies 
carried out by IMarschall et alj ([19881 11989ft confirmed 
IM Vir to be both double-lined and eclipsing, and to be 
composed of a G7 star and a late-K or early-M star, thus 
piquing our interest. The orbital period was estimated 
as 1.3085 days. 

Very little data on this binary have been published 
since, other than sparse photometry and occasional re- 
ports on the chromosph e ric activity and X-ray flaring 
([Strassmeier et~aT1 Il993l: IPandev fc Shier] [2001 . We 
present here extensive spectroscopic and differential pho- 
tometric measurements that allow us to determine its 
fundamental properties with very high accuracy. The 
observations are presented in §[5]and §[31 followed by a de- 
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tailed light-curve analysis in §U accounting for the pres- 
ence of spots. The absolute dimensions are derived in 
§[5] with careful consideration of potential sources of sys- 
tematic error, essential for the results to be useful. The 
comparison with models of stellar evolution is discussed 
in §|51 and we summarize our conclusions in §[7] 

2. SPECTROSCOPIC OBSERVATIONS AND REDUCTIONS 

IM Vir was observed at the Harvard-Smithsonian Cen- 
ter for Astrophysics (CfA) with three nearly identical 
echelle spectrographs on the 1.5-m Tillinghast reflector 
at the F. L. Whipple Observatory (Mount Hopkins, Ari- 
zona), the 1.5-m Wyeth reflector at the Oak Ridge Obser- 
vatory (Harvard, Massachusetts), and the 4.5-m equiva- 
lent Multiple Mirror Telescope (also on Mount Hopkins, 
Arizona) prior to its conversion to a 6.5-m monolithic 
telescope. Photon-count ing intensified Retico n detectors 
('Digital Speedometers': iLathaml Il985l 11992T ) were used 
in each case to record a single 45 A echelle order cen- 
tered at a wavelength of 5188.5 A, featuring the gravity- 
sensitive lines of the Mg I b triplet. The resolving power 
provided by this setup is A/AA » 35,000. Nominal 
signal-to-noise (S/N) ratios for the 138 spectra we ob- 
tained range from 13 to 58 per resolution element of 
8.5 kms -1 . The first observation was taken in 1984 Jan- 
uary 1, and monitoring continued until 2009 May 10. A 
handful of the e arly spectra are t he same ones included 
in the work of iSilva et all (| 19871 ) , who discovered the 
radial-velocity variability, but have been re-reduced and 
analyzed here with much improved methods, as we now 
describe. 

All of our spectra are double-lined, but the secondary 
is comparatively quite faint, accounting for only 6% 
of the light of the primary (see below). Radial ve- 
locities for both stars were derived using TODCOR, a 
tw o-dimensional cross-corr elation technique introduced 
bv lZucker fc Mazehl (|1994l ). This method uses two tem- 
plates, one for each component of the binary, which we 
selected from a large library of synthetic sp ectra based on 
mode l atmospheres by R. L. Kurucz (see lLatham et all 
2002). These templates have been calculated for a wide 
range of effective temperatures (T e g), surface gravities 
(log<?), rotational velocities (vsmi when se en in projec- 
tion), and metallicities ([m/H]). Following iTorres et alJ 
(2002), the optimum templates were selected by run- 
ning extensive grids of cross-correlations with TODCOR, 
seeking to maximize the average correlation weighted by 
the strength of each exposure. Because of the faint- 
ness of the secondary, only the parameters for the bright 
star can be determined reliably from our spectra. Ini- 
tially we assumed logg = 4.5, and determined T e g and 
■usinz for a range of metallicities between [m/H] = —1.0 
and [m/H] = +0.5. The optimal values were found 
by interpolation in [m/H]. We then repeated the cor- 
relations for a lower value of logg = 4.0 in order to 
bracket the estimate from the analysis described later 
(log<7 = 4.379). Interpolation in log<? to this final value 
resulted in T cS = 5570±100K and wsin i = 43±2kms -1 . 
The formal mctallicity that maximizes the correlation is 
[m/H] = —0.1. While this suggests an overall abundance 
close to solar, the uncertainty is likely to be at least 
0.25 dex, and we consider the estimate primarily as a 
free parameter included to optimize the match between 



the synthetic templates and the observed spectra. We 
note, furthermore, that because of our narrow spectral 
window, there is a strong correlation between tempera- 
ture and metallicity. In this case, however, the T e ff value 
reported above is consistent with various photometric es- 
timates described below in §[5j in turn lending more cre- 
dence to the metallicity estimate. For the final velocity 
determinations we adopted the set of template parame- 
ters closest to the values above that maximizes the av- 
erage correlation (T e g = 5750 K, vsmi — 40 kms -1 , 
log<? = 4.5, and solar composition). Small differences 
in these template parameters compared to the interpo- 
lated values do not affect the velocities significantly. For 
the secondary template we adopted the same metallic- 
ity, and parameters consistent with the results from §[5j 
T off = 4250 K, vsmi = 25 kms" 1 , logg = 4.5. 

In addition to the velocities, we have determined the 
light ratio between the components at the mean wave- 
length of our observations, following iZucker fe Mazehl 
(1994), accounting for the difference in line blocking be- 
tween the primary and the much cooler secondary. 5 The 
value we obtain is 1^1 1\ = 0.06 ± 0.01. 

Although TODCOR significantly reduces systematic 
errors in the radial velocities caused by line blending, 
residual effects can remain as a result of shifts of the 
spectral lines in and out of our narrow spectral window 
as a function of orbital phase. Experience has shown 
that these effects must be examined on a case-by-case 
basis. We investigated them here by mean s of numerical 
simula tions sim ilar to those descr i bed b y lLatham et alJ 
(p96Tl (see also ITorres etalMMA f2000h . We generated 
synthetic composite spectra matching our observations 
by combining copies of the primary and secondary tem- 
plates used above, shifted to the appropriate velocities 
for each of the exposures as predicted by a preliminary 
orbital solution, and scaled to the observed light ratio. 
These synthetic observations were then processed with 
TODCOR in exactly the same way as the real spectra, 
and the resulting velocities were compared with the input 
shifts. The differences for IM Vir are shown in Figure [TJ 
and are less than 0.5 kms -1 for the primary but reach 
values as large as 13 kms -1 for the secondary. Similarly 
large differences have been found occasionally for other 
systems using the same instrumentation (e.g., AD Boo , 
GX Gem, VZ Cep:lclaus en et alJl2?M ILacv et al.H 200a 
ITorres fc Lacvll2009l) . We have applied these differences 
as corrections to the raw velocities of IM Vir. The impact 
on the masses, however, is very small: only 0.26% for the 
primary and 0.14% for the secondary. The reason for this 
is that the large corrections for the secondary are simi- 
lar and of the same sign at both quadratures, therefore 
amounting mostly to an overall systematic shift rather 
than a change in the velocity semi-amplitude. Similar 
adjustments based on the same simulations were made 
to the light ratio, and are already included in the value 
reported above. 

The stability of the zero-point of the CfA velocity sys- 
tem was monitored by means of exposures of the dusk 
and dawn sky, and small run-to-ru n correct i ons w ere ap- 
plied in the manner described by |Lathairj 11992J). The 

5 This is to correct the ratio of continuum heights for the fact 
that the lines of the secondary are intrinsically stronger, and sub- 
tract proportionately more flux from the spectral window. 



IM Vir 



3 



TABLE 1 

Radial velocity measurements of IM Vm. 



HJD 


RVi 


RV 2 


01 


C2 


(O-Ch 


(o-c) 2 




(2,400,000+) 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 


(kms 1 ) 


(kms -1 ) 


Phase 


45719.9949 .... 


-67.58 


+144.06 


1.54 


13.01 


+0.30 


+9.61 


0.1661 


45743.0126 .... 


+105.93 


-116.81 


0.83 


6.98 


+1.08 


+3.76 


0.7555 


45750.8945 .... 


+100.91 


-122.28 


1.58 


13.30 


-2.51 


-3.82 


0.7785 


45754.0581 .... 


-74.09 


+148.80 


1.54 


12.97 


+1.10 


+3.56 


0.1961 


45754.0603 .... 


-71.28 


+149.09 


1.90 


16.05 


+4.23 


+3.38 


0.1977 



Note. — This table is available in its entirety in 
A portion is shown here for guidance regarding its 



a machine-readable form in the online journal, 
form and content. 
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Fig. 1. — Systematic corrections applied to the raw radial veloc- 
ity measurements as a function of orbital phase (top) and radial 
velocity (bottom). Filled symbols represent the primary, and open 
symbols are for the secondary. 

final velocities, including these offsets as well as the cor- 
rections for systematics, are listed in Table [T] along with 
their uncertainties. The median values for the velocity 
errors are 0.90 kms -1 for the primary and 7.6 kms -1 for 
the secondary. 

Preliminary single-lined orbital solutions carried out 
separately for the primary and secondary indicated a sig- 
nificant zero-point difference (i.e., a difference in the sys- 
temic velocity 7), which is often seen by many investiga- 
tors in cases where there is a slight mismatch between the 
spectra of the real stars and the template s used for the 
cross -correlations (see, e.g., [Popper 2000; Griffi n~et al.1 
2000). Numerous experiments with other templates did 
not remove the offset, which is about 4 kms -1 (the sec- 
ondary velocities being systematically larger). In IM Vir 
this most likely arises because of the cool temperature of 
the secondary, and the fact that synthetic templates for 
such stars become increasingly unrealistic due to miss- 
ing opacity sources in the models. In order to prevent 
this offset from affecting the velocity semiamplitudes, we 
have therefore accounted for it by including it as an ad- 
ditional free parameter in our double-lined solution. Ef- 
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Fig. 2. — Phase-folded radial velocity observations for IM Vir 
(filled circles for the primary, open circles for the secondary), along 
with our best-fit model. Phase 0.0 corresponds to primary eclipse. 
The residuals from the fit are shown in the bottom panels. 

fectively this means we allow each component to have its 
own center-of-mass velocity, without significantly affect- 
ing the velocity semiamplitude of the component, and 
hence without affecting the individual mass determina- 
tion. On the other hand, the offset implies additional un- 
certainty in the true center-of-mass velocity of the binary, 
beyond the formal errors listed below. We do not expect 
a significant template mismatch for the solar-type pri- 
mary, so the systemic velocity is most likely to be closer 
to the value for that star. 

The results are presented in Tabled in which the mea- 
surements have been weighted according to their uncer- 
tainties in the usual way, and the errors rescaled by itera- 
tions to achieve a reduced \ 2 value near unity, separately 
for the primary and secondary. No indication of signifi- 
cant eccentricity was found, as expected for such a short 
orbital period, and the final orbit was therefore consid- 
ered circular. The observations and computed curve are 
displayed graphically in Figure^ along with the residu- 
als, which are also listed in Table [TJ 

3. PHOTOMETRIC OBSERVATIONS 
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TABLE 2 

Spectroscopic orbital solution for IM Vir. 

Parameter Value 

Adjusted quantities 

P (days) 1.30861497 ± 0.00000034 

Ti (HJD-2,400,000) a 52402.87420 ± 0.00052 

K\ (kms" 1 ) 92.684 ± 0.088 

K 2 (kms" 1 ) 136.84 ± 0.74 

e b 

7 (kms" 1 ) +12.221 ± 0.082 

ARV (kms _1 ) d -3.96 ± 0.69 

Derived quantities 

Mi sin 3 i (M Q ) 0.978 ± 0.012 

M 2 sin 3 i (M Q ) 0.6621 ± 0.0044 

q = Ms /Mi 0.6773 ± 0.0037 

ttl sin i (10 6 km) 1.6678 ± 0.0016 

a 2 sin i (10 6 km) 2.462 ± 0.013 

asini (10 6 km) 4.130 ± 0.013 

Other quantities pertaining to the fit 

<Jl (kms" 1 ) 0.99 

CT 2 (kms -1 ) 8.39 

AUs 138 

Time span (days) 9241.7 

a Time of primary eclipse. 
b Fixed; see text. 

c True uncertainty may be larger due to systematic effects; sec 
text. 

d Zero-point difference between the primary and secondary ve- 
locities, in the sense primary minus secondary. 

Differential photometric observations of IM Vir were 
conducted in 2006 April and May using the 0.4-m 
Ealing Cassegrain reflector of the Gettysburg College 
Observatory (GCO, Gettysburg, Pennsylvania). The 
camera was a Photometries (Roper Scientific) CH-350 
thermoelectrically-cooled unit equipped with a SITc 
003B back-illuminated scientific grade CCD chip and 
Bessell BVRI filters. Exposure times were 120, 60, 25, 
and 20 seconds, respectively. We obtained 597, 672, 630, 
and 631 observations in B, V, R, and /, with a ca- 
dence of about 6 minutes. The field of view of the GCO 
camera is approximately 18 arcmin, and since IM Vir 
is a relatively bright target, this severely limited our 
choice of reference stars of comparable magnitude and 
color to two: BD-05 3573 ('comp', a = 12 h 48 m 36?40, 
5 = -5°53'33."6, J2000.0; V = 10.39, B-V = 0.67) and 
HD 111427 ('check', a = 12 h 49 m 14?94, S = -5°49'20."7, 
J2000.0; V = 9.40, B-V = 0.67). The colors of these 
two stars are in fact nearly identical to that of the vari- 
able, which is B — V = 0.66 (see §0. 

Differential photometry was performed on IM Vir and 
the two reference star s in all our images by mean s 
of MIRA-AP software (http:/ /www.mirametrics.com/|). 
We employed standard aperture photometry techniques 
to derive instrumental magnitudes, setting the radius of 
the measuring apertures for each night using a standard 
value of 2.5 times the FWHM of the seeing disk, based 
on previous curve-of-growth calibrations using the same 
equipment. Typical errors as represented by the scat- 
ter of the comp— check differences are 0.0132 mag in B, 
0.0124 mag in V, 0.0135 mag in R, and 0.0150 mag in /. 
The primary eclipse is ^0.62 mag deep in V, while the 
depth of the secondary is only ~0.08 mag. 

Examination of the raw photometry revealed slight 
trends in the comp— check differences that changed from 
night to night and are most likely of instrumental origin. 



TABLE 3 

Differential B-band photometry of IM Vir. 



HJD 

(2,400,000+) 


Phase 


Comp-Var 
(mag) 


Comp-Check 
(mag) 


53835.58358 ... 


0.8288 


0.811 


0.966 


53835.58732 . . . 


0.8317 


0.849 


0.995 


53835.59108 . . . 


0.8345 


0.852 


1.000 


53835.59482 ... 


0.8374 


0.853 


1.007 


53835.59857 . . . 


0.8403 


0.854 


1.005 



Note. — This table is available in its entirety in a 
machine-readable form in the online journal. A portion is 
shown here for guidance regarding its form and content. 



TABLE 4 

Differential V-band photometry of IM Vir. 



HJD 

(2,400,000+) 


Phase 


Comp-Var 
(mag) 


Comp-Check 
(mag) 


53835.58481 ... 


0.8298 


0.848 


0.994 


53835.58855 . . . 


0.8326 


0.846 


0.999 


53835.59230 . . . 


0.8355 


0.839 


0.989 


53835.59605 . . . 


0.8383 


0.838 


0.998 


53835.60705 . . . 


0.8468 


0.826 


0.985 



Note. — This table is available in its entirety in a 
machine-readable form in the online journal. A portion is 
shown here for guidance regarding its form and content. 



TABLE 5 

Differential R-band photometry of IM Vir. 



HJD 

(2,400,000+) 


Phase 


Comp-Var 
(mag) 


Comp-Check 
(mag) 


53835.62654 . . . 


0.8616 


0.584 


1.027 


53835.63031 . . . 


0.8645 


0.660 


1.011 


53835.63406 . . . 


0.8674 


0.685 


1.002 


53835.63781 . . . 


0.8703 


0.704 


1.003 


53835.64157 . . . 


0.8731 


0.708 


1.005 



Note. — This table is available in its entirety in a 
machine-readable form in the online journal. A portion is 
shown here for guidance regarding its form and content. 



TABLE 6 

Differential /-band photometry of IM Vir. 



HJD 

(2,400,000+) 


Phase 


Comp-Var 
(mag) 


Comp-Check 
(mag) 


53835.58592 ... 


0.8306 


0.933 


1.013 


53835.58967 . . . 


0.8335 


0.933 


1.014 


53835.59342 . . . 


0.8363 


0.939 


1.034 


53835.59717 ... 


0.8392 


0.929 


1.025 


53835.60091 . . . 


0.8421 


0.946 


1.032 



Note. — This table is available in its entirety in a 
machine-readable form in the online journal. A portion is 
shown here for guidance regarding its form and content. 



Neither star is known to be variable, and we find no peri- 
odicities or other discernible patterns. Consequently, we 
de-trended the IM Vir measurements by removing the 
median value of the comp— check differences calculated 
over intervals of a few hours. Even after this correction, 
the sparser second half of the data presents a system- 
atic ~0.01 mag difference (fainter) in the out-of-eclipse 
light level compared to the first half, and other occasional 



IM Vir 



5 



TABLE 7 
Times of eclipse for IM Vir. 



HJD 








O-C 




(2,400,000+) 


(days) 


Type a 


Instr. 


(days) 


Source 


46042.99984 


0.00097 


I 


ccd 


-0.00561 


1 


51274.8750 




1 


ccd 


+0.02690 


2 


51885.9610 




1 


vis 


-0.01029 


2 


52402. 87420 b 


0.00052 


1 


spec 


0.0 


4 


53474.62948 


0.00055 


I 


ccd 


-0.00038 


3 


53843.65995 


0.00010 


1 


ccd 


+0.00067 


4 


53845.62427 


0.00052 


11 


ccd 


+0.00207 


1 


53877.68325 


0.00038 


I 


ccd 


-0.00002 


4 



a Eclipses labeled as T for primary, and 'IF for secondary. 

b Mean epoch from the radial vel ocities. 

c Source: 1. Mean epoch from the Manfroid et al. (1991) photom- 
etry; 2. B.R.N.O. database | |http://var. astro. cz\ ; 3. IQgloza et alj 
(2008); 4. This work. 



changes possibly due to temporal evolution of surface in- 
homogeneities (spots) on the surface of one or both stars. 
This would indicate a timescale for evolution of these fea- 
tures of a few weeks. For these reasons we have chosen to 
restrict our analysis to the first half of the observations 
(spanning 22 days), which provide complete coverage of 
both eclipses and are cleaner overall. The complete data 
set is given for all filters in Tables [3HS] in its original form, 
i.e., without the nightly corrections, and we list also the 
comp— check differences. 

Additional photometry of IM Vir has been reported 
by iManfroid et al.1 (|199lD in the Stromgren uvby sys- 
tem. These data were obtained some twenty years earlier 
than our own measurements, between 1983 May and 1986 
July. Unfortunately the coverage of the eclipses is very 
incomplete, as shown later, so these data are not useful 
for determining the geometric parameters of the system. 
Nevertheless, it is possible to extract an average time of 
eclipse, which we present below, as well as brightness ra- 
tios in the different bands, which will be used in §[5] to 
deconvolve the light of the two stars in order to obtain 
photometric estimates of the temperatures and metallic- 
ity 

3.1. Ephemeris and times of minimum 

The ephemeris adopted for the remainder of this paper 
is the one calculated from the spectroscopy, 

Min I = 2,452,402.87420(52) + 1.30861497(34).E , (1) 

which, by virtue of the 25-year radial-velocity coverage, 
is much more accurate than could ever be obtained from 
the 64-day interval of the photometric observations. The 
uncertainties are indicated above in parentheses. 

Three times of eclipse from the BVRI photometry have 
been measured by fitting the light-curve model described 
below (including spot parameters) to each night with suf- 
ficient coverage of a primary or secondary minimum, si- 
multaneously in all four passbands. The only adjustable 
parameter allowed in these fits wa s a time shift. A s imilar 
procedure was followed with the IManfroid et~aTl (|1991f ) 
photometry, with the primary luminosity in each band 
added as a free parameter and no spots considered. The 
resulting eclipse timings for IM Vir are collected in Ta- 
ble [71 along with an average time of eclipse from the 
spectroscopy, and the few additional times found in the 
literature. 



4. LIGHT CURVE ANALYSIS 

The differential photometry described in §[3]was phase- 
folded with the adopted ephemeris, and analyzed with 
the 2003 version o f the WD code (|Wilson fe Devinnevl 
[19711: IWilsoni ri979). For the reasons indicated earlier, 
only the first half of the data were used (22 nights worth). 
The main light curve parameters adjusted are the or- 
bital inclination (i), the pseudo-potentials (f^i and f^), 
the secondary effective temperature (T^ff^), the primary 
luminosity (Li), and a phase shift. The latter would 
normally be unnecessary given the accuracy of the spec- 
troscopic ephemeris and the short interval of the photo- 
metric observations. However, as we describe below, the 
light curves of IM Vir present distortions we ascribe to 
spots, which can cause asymmetries in the eclipses that 
may appear as phase shifts. Since these distortions are 
modeled here only in a simplified way, it is prudent to 
include a phase shift as an additional parameter to sup- 
plement the spot model. The primary temperature was 
held fixed at the spectroscopic value of 5570 K, since the 
light curves do not provide a strong enough constraint 
on both temperatures, but only on their ratio. Emergent 
intensities used in the pr ogram were taken from mode l 
atmospheres described by Ivan Hamme fe Wilson! |2003). 
Square-root limb-darkening coefficients for the Johnson- 
Co usins BVRI ba nds were interpolated from the tables 
bv iClaretl (|2000a| ). and adjusted dynamically according 
to the current temperatures and surface gravities of the 
stars at each iteration. The reflection albedos were fixed 
at the value 0.5, appropriate for stars with convective en- 
velopes, and the gravity darkening exponents were set to 
0.3 4 for the prima ry and 0.40 for the secondary, follow- 
ing [claret] (|2000bf ). The orbit was initially considered to 
be circular, and spin-orbit synchronization was assumed 
based on the evidence presented in §[5] The photometric 
measurements in the four passbands were fitted simulta- 
neously. Outliers were rejected by 3-cr clipping based on 
preliminary solutions. Convergence in the final fit was 
considered to have been achieved when the corrections 
to the elements were smaller than the internal errors in 
three consecutive iterations. 

As is common in eclipsing binaries with late-type com- 
ponents, the light curves of IM Vir show out-of-eclipse 
modulations that are attributable to spots. There is also 
a slight difference in the light level at the two quadra- 
tures. These variations can in principle be modeled with 
the WD code in the approximation of circular spots of 
uniform temperature. In practice, however, such mod- 
eling is fraught with difficulties and there is abundant 
literature discussing problems of indeterminacy and non- 
uniqueness, p articularly when u sing data of limited qual- 
ity (see, e.g.. lEkerl I1996L [19991 and references therein). 
Nevertheless, because we are interested in recovering the 
geometric parameters of the system as free as possible 
from systematic errors, we have made an effort to ac- 
count for these modulations by considering up to two 
spots. 

While the subset of the photometric observations used 
here spans only 22 days, our spectroscopic coverage ex- 
tends for more than 25 years. Given the changes in the 
light curves mentioned in §[3j the photometry is modeled 
separately from the radial velocities as any spot model 
would not apply to both. In a detached system such as 
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TABLE 8 

Light curve parameters of IM Vir from our simultaneous 
BVRI fit. 

Parameter Primary Secondary 

Geometric properties 

Phase shift 0.0006 ± 0.0001 

i (deg) 87.24 ± 0.16 

Q 6.298 ± 0.023 7.081 ± 0.038 

r pole 0.1776 ± 0.0007 0.1145 ± 0.0007 

r point 0.1796 ± 0.0008 0.1151 ± 0.0008 

r alde 0.1784 ± 0.0007 0.1147 ± 0.0007 

r back 0.1776 ± 0.0008 0.1150 ± 0.0007 

r vol a 0.1785 ± 0.0008 0.1146 ± 0.0010 

Radiative properties 

T cff (K) 5570 b 4246 ± 16 

L 2 /ii B band 0.04805 ± 0.00008 

L2/L1 V band 0.07499 ± 0.00016 

L2/L1 R band 0.10960 ± 0.00027 

L2/L1 / band 0.14267 ± 0.00036 

Albedo c 0.5 

Gravity darkening .... 0.34 0.40 

Limb darkening coefficients (square root law) 

x B band 0.625 0.904 

y B band 0.240 -0.057 

x V band 0.364 0.644 

y V band 0.450 0.190 

x R band 0.237 0.462 

y R band 0.517 0.354 

x I band 0.137 0.269 

y I band 0.539 0.475 

Spots properties 

Latitude (deg) -60 -60 

Longitude (deg) 307.3 ± 4.9 331.8 ± 4.7 

Radius (deg) 26.4 ± 1.2 36.0 ± 2.5 

T cff factor c 0.95 0.80 

Residuals and number of observations 

<7g / N B 0.01347 / 438 

a v / N v 0.01295 / 495 

o R I N a 0.01300 / 443 

07 / Nj 0.01399 / 455 

a Volume radius derived from the fitted parameters. 

Fixed according to the spectroscopic analysis. 
c Fixed; see text. 

this the photometry does not constrain the mass ratio, 
so the value of q = 0.6773 was adopted from the spec- 
troscopy. 

The parameters representing the spots in the WD 
model are the latitude (6* s ), longitude (</> s ), angular radius 
(r s ) and temperature contrast relative to the photosphere 
(T e ff jS /T e ff). Only cool spots have been considered. The 
latitude is essentially unconstrained by our data, and the 
spot size and temperature factor are strongly correlated 
with each other and cannot usually be adjusted simulta- 
neously. Thus, only the longitude and the spot size were 
adjusted at the same time as the geometric parameters, 
while the latitude and contrast factor were held fixed and 
estimated through a grid search. We explored values over 
a wide range in 8 S from +80° to —80° in steps of 20°, a 
range in the primary T e g iS /T c f[ from 0.85 to 0.95, and a 
range from 0.75 to 0.85 for the secondary T e s yS /T e s, both 
with a step size of 0.05. As indicated above, it is virtu- 
ally impossible to tell which star has the spot(s), so we 
have considered here only three possible scenarios: two 
spots on the primary, two on the secondary, or one spot 
on each star. 

The best solutions were typically obtained with the 
spots located in the southern hemisphere, although we 
do not assign any particular physical significance to this 
as we only consider the spot modeling as a means of re- 
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Fig. 3. — BVRI observations of IM Vir along with our best-fit 
model. The curves are shifted vertically for clarity. Photometric 
residuals are shown in the bottom panels, in the same order as the 
top curves. 
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Fig. 4. — Enlargement of the eclipse phases in Figure [3] Note 
the different vertical scales for the primary and secondary eclipses. 

moving a perturbation from the underlying light curve. 
The overall best fit has one spot on each component, and 
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TABLE 9 

Summary of single-passband light curve fits, and estimated uncertainties. 



Individual light-curve fits 



Standard error estimates 



Parameter 



B 



V 



ft 



""Bands ""Spots 



"WD 



Adopted fit 
Value cr 



Phase shift . 

i (deg) 

Teff,2 

Qi 

n 2 

r vo l,i 

r vo l,2 ■ 
L 2 /Li 
L2/L1 

L 2 /L 1 



(B) 
(V) 
(R). 
{!)■ 



0.0007 
87.38 
4272 
6.257 
7.091 
0.1798 
0.1144 
0.04999 



0.0007 
87.22 
4191 
6.289 
7.080 
0.1788 
0.1146 

0.05366 



0.0006 
86.91 
4231 
6.304 
7.003 
0.1783 
0.1161 



0.10928 



0.0004 
87.83 
4304 
6.431 
7.220 
0.1743 
0.1119 



0.15179 



0.00014 
0.38 
49 
0.077 
0.090 
0.0024 
0.0017 



0.00010 
0.08 
56 
0.018 
0.044 
0.0006 
0.0009 
0.0041 
0.0054 
0.0065 
0.0068 



0.00010 
0.16 
16 
0.023 
0.038 
0.0008 
0.0010 
0.00021 
0.00030 
0.00042 
0.00050 



0.0006 
87.24 
4246 
6.298 
7.08 
0.1785 
0.1146 
0.0480 
0.0750 
0.1096 
0.1427 



0.0002 
0.42 
83 
0.082 
0.11 
0.0026 
0.0022 
0.0041 
0.0054 
0.0065 
0.0068 



g 

< 0.04 
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Fig. 5. — Effect of the spots on the light curves. The variations 
shown correspond to the adopted scenario with one cool spot on 
each star. 

a reduced x 2 square that is only marginally better than 
the other two scenarios (1.5% lower than the case with 
two spots on the primary, and 3.4% lower than the two- 
spot configuration on the secondary). The parameters 
of this best fit are listed in Table [Sj and the synthetic 
curves are shown together with the observations in Fig- 
ure [3] and Figure 2] Among the quantities derived from 
this fit, the light ratios allow for an important consis- 
tency check against the ratio obtained directly from our 
spectra. Interpolating between B and V to the mean 
wavelength of our spectroscopic observations, we obtain 
^2/^1 = 0.066 ± 0.005, which agrees well with the value 
of £ 2 /£i = 0.06 ±0.01 from §[2 The effect of the spots on 
the light curves is illustrated in Figure [51 and is seen to 
range from ^0.025 mag in B to ^0.015 mag in / (peak- 
to-peak). A depiction of the location of the spots on each 
star is shown in Figured! The spot on the primary cov- 
ers 5.2% of its surface, while that of the secondary covers 
9.5%. 

In order to provide more realistic uncertainties for the 
geometric and radiative parameters than the internal er- 
rors reported in Table we have considered additional 
sources of error as follows: (a) We ran WD solutions sep- 
arately for each passband, with the spot parameters fixed 
at their final values, and examined the differences in the 
parameters. These fits are summarized in Table |H1 The 
dispersion in each parameter about the average of the 
four bands (ffBands) was taken as an additional contri- 
bution to the overall uncertainty, (b) We calculated the 
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Fig. 6. — Spot location on each star as viewed from the Earth 
at different orbital phases, in the adopted scenario in which each 
component has one cool spot. The stars and their separation are 
rendered to scale. 



range in each parameter from the simultaneous BVRI 
fits in the three spot scenarios, and adopted half of this 
range as an additional contribution to account for the 
degeneracy in the spot modeling, (c) As a check on the 
internal errors from WD, we continued the iterations in 
our adopted fit beyond convergence for another 200 steps, 
and we examined the scatter of those 200 solutions. For 
the main parameters (i, T c g^, ^1, ^2) the scatter was 
only a small fraction of the internal errors, but for the 
light ratios it was typically a factor of two larger. We 
adopted the larger of the two estimates in each case. 

These three sources of error (cBands, cspots) and ctwd 
in Tabled were combined quadratically, and are the ones 
we assign to the final light-curve parameters. Changes 
in the adopted mass ratio within its error have no addi- 
tional effect. The parameters and these errors, which we 
believe to be realistic, are listed in the last two columns 
of Table [9] Further solutions were carried out allowing 
the eccentricity to vary, but in all cases we found the 
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Fig. 7. — Differential BVRI photometry of IM Vir corresponding 
to the second half of our data set, which we have not included in 
the analysis for the reasons described in the text. The solid curves 
are the same best-fit model shown in Figure[3] (without the spots), 
adjusted for a slight overall brightness change and a phase shift. 

result to be insignificant compared to its error, consis- 
tent with the indications from spectroscopy. Third light 
was also tested for at various stages of the analysis, but 
was always found to either converge toward negative (un- 
physical) values, or to be consistent with zero. Different 
limb-darkening laws were tested as well (linear, logarith- 
mic), but the differences with the results in Table [5] were 
minimal (well within the errors). 

The final solution indicates the two stars are nearly 
spherical, the difference between r po i nt and r po i c being 
only 1.1% for the primary and 0.5% for the secondary. 
The primary eclipse is annular (48% of the light of that 
star blocked), while the secondary is total (see Figure |4|), 
with the totality phase lasting 28 minutes. 

The second half of our photometry, not used here be- 
cause of likely changes in the spots and more incomplete 
phase coverage, is compared with our final model for the 
geometry in Figure excluding the spot terms. The 
larger scatter is obvious, but the fit is still quite reason- 
able. Figure [8] displays our fits to the sparser Stromgren 
photometry bv lManfroid et al. (1991). Brightness ratios 
from these solutions (0.0167, 0.0353, 0.0599, and 0.0775 
in uvby, respectively) are used below to separate the light 
of the two components and derive a photometric estimate 
of their effective temperatures and metallicity. 

5. ABSOLUTE DIMENSIONS 

The absolute masses and radii for the components of 
IM Vir follow from our spectroscopic and light curve solu- 
tions described in the preceding sections. The effect that 
spots may have on the relative radii has been consid- 
ered explicitly in our light-curve modeling, and included 
in our error estimates. Similar effects may influence the 
masses, although we do not expect them to be signifi- 
cant since they would tend to average out over the ex- 
tended coverage of our spectroscopy. Nevertheless, the 
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Fig. 8. — Stromgren uvby photo metry on the standard system 
published by Manfroid et al. (1991), but not used in our analysis, 
compared with our best-fit model. 

spot models in §[4] allow us to quantify the importance of 
these effects. Figure |9] shows the distortions in the radial 
velocities expected from each of the three spot scenarios 
considered earlier. As a test, we perturbed the radial 
velocities using these curves and repeated the spectro- 
scopic orbital solutions. Compared to our adopted fit 
(Tabled, the differences were less than 0.6% in the min- 
imum masses. The differences in asinz, which affects 
the absolute radii, were less than 0.2%. To be conser- 
vative, for the calculation of the final mass and radius 
uncertainties we have augmented the errors in Table [2] 
by adding in quadrature half of the range in M12 sin 3 i 
and asini obtained from the three spot configurations. 
With this, the absolute masses of IM Vir are determined 
to 1.2% and 0.7% for the primary and secondary, and the 
absolute radii to 1.5% and 1.9%. 

Next in importance to the masses and radii are the 
effective temperatures of the stars. The temperature 
difference (or ratio) is accurately determined from the 
light curves, but the absolute scale is set by the pri- 
mary value, which is fixed in our analysis. As a check 
on the spectroscopic value of 5570 K adopted here, we 
made use of absolute photometry available for IM Vir 
from a variety of sources, along with color/temperature 
calibrations. Brightness measurements we re compiled 
from the 2 MASS catalog jlCutri et al.ll2003f) . the Tycho- 
2 catalog (lH0g et a l. 2000), Stromgren photom etry from 
iManfroid et all (|199lb and lMorale et all (|1996h . and our 
own BV measurements outside of eclipse, transformed 
to the standa rd system using the co mparison and check 
stars. For the Ma nfroid et all (|199lD photometry we con- 
sidered only t he da ta out of eclipse. The Tycho-2 and 
iMorale et alJ (j 19961 ) measurements are quite numerous 
and can be considered to be little affected by the min- 
ima. The date of the 2MASS measurements indicates 
they were taken during the egress of a secondary eclipse; 
corrections of order 0.03 mag (slightly different in each 
band) were applied to JHK S based on our light curve 
solutions, extended to the near-infrared. The adopted 
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TABLE 10 

OUT-OF-ECLIPSE COMBINED-LIGHT 
MAGNITUDES AND COLORS OF IM VlR. 



Passband 


Value 


Magnitudes 




JD 


10.234 ± U.030 


V 


y.0/4 ± U.030 


L 






y 


O *7A O X (1 (117 

y. i 4o ± U.U1 ( 




in /iqq _i_ n nQO 

iu.4oo ± u.uoy 


V T 


y.768 ± 0.030 


J 


8.176 ± 0.020 a 


H 


7.712 ± 0.025 a 


K a 


7.634 ± 0.024 a 


Color indices 




B — V 


0.661 ± 0.028 


b-y 


0.465 ± 0.032 


V — J 


1.372 ± 0.028 


V-H 


1.832 ± 0.030 


V - K s 


1.905 ± 0.031 


B T - V T 


0.715 ± 0.057 


Vt-Ks 


2.099 ± 0.038 


J — K s 


0.533 ± 0.031 



a Corrections of 0.026, 0.030, and 
0.035 mag have been applied to the mea- 
sured J HK S values to account for the fact 
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Fig. 9. — Effect of the spots on the radial velocities, for the 
three scenarios explored with our WD solutions (solid lines for the 
primary, dashed for the secondary). 

magnitudes for the combined light in all passbands are 
listed in Table [TOl With these we formed 8 different color 
indices, also listed in the table, and then deconvolved 
the light to obtain the individual indices for the two 
stars. The light ratios in B and V for the deconvolution 
were adopted from Table [H and those for the Stromgren 
passbands from our fits in the previous se ction. For the 
other passbands we used model isoch rones (iGirardi et al" 
120001 for the Tycho-2 bands, and IBaraffe et alj I199S 
for the JHK S bands, after conversion to the 2MASS 
system). Two different isochrone ages were tried 
(1 Gyr and 5 Gyr), although the differences were well 
within the errors. Co l or /te m perature calibratio n s from 
iRaimrez fe Melendezl (120051). ICasagrande et ail (|2006l ). 
and iGonzalez Hernandez fc Bonifacio! (|2009f ) were then 
applied to the primary indices, assuming solar metal- 
licity. Although the color indices are not independent 
of each other, they do serve to test the internal consis- 
tency between the various photometric systems and the 
three different calibrations. The interagreement is very 
good (typically within 150 K). The weighted average is 
5560 ± 100 K when using the 1 Gyr isochrones for the 
light ratios, and 40 K less when adopting 5 Gyr. We 
regard these to be in very good agreement with the spec- 
troscopic value, considering the uncertainties. Changing 
the metallicity for the calibrations to the rough estimate 
from §[5] ([m/H] w —0.1) reduces the photometric tem- 
peratures by less than 10 K. A similar exercise using 
the indices for the sec ondary and an add i tional calibra- 
tion for cool stars by ICasagrande et al.l (|2008f ) gives a 
more uncertain result of 4380 ± 220 K (as expected from 
the faintness of the star), but still consistent with the 
much more accurate value based on the light curve anal- 



obscrvation was 
during the egress of a secondary eclipse 
(sec text). 

ysis, which is 4250 ± 130 K. These tests suggest that the 
adopted temperatures for IM Vir are accurate. They cor- 
respond to spectral types of G7 V for the primary and 
K7 V for the secondary. 

An additional quantity of great importance for the in- 
terpretation of the masses, radii, and temperatures is 
the chemical composition. Beyond our estimate in §[2] of 
[m/H] w —0.10 ± 0.25, a spectroscopic estim ate from a 
compo site spectrum of IM Vir was reported bv lDall et alJ 
(|2007f) as [Fe/H] = -0.53 ± 0.16, based on an effective 
temperature some 200 K cooler and a higher surface 
gravity than we derive for the primary, which dominates 
the light. It is unclear how accurate this determination 
is, in view of those differences. A photometric estimate 
for the primary may be derived from t he out-of-eclipse 
Strom gren measurements of IM Vir by IManfroid et al.l 
(|1991l) and lMorale et all (|1996l ). after removing the light 
contribution of the secondary using the light ratios ob- 
tained in §|U The metallicity relation by Holmbcrg et alJ 
(2007) gives a rather poorly determined value of [Fc/H] 
~ —0.37 ±0.47, in which the uncertainty includes photo- 
metric errors as well as the scatter of the calibration. The 
secondary is too cool for this calibration and other similar 
ones based on Stromgren in dices, but is w i thin r ange of 
the near- infrared formula bv lBonfils et al.l (|2005[ ). which 
yields [Fe/H] w -0.26 ± 0.26. This again includes all 
photometric error s and the scatter of the calibration. It 
has been noted by Uohnson fc Appsl (|2009D that this lat- 
ter color/temperature relation appears to underestimate 
the metallicity of late-type stars of solar composition or 
greater, by approximately 0.3 dex; for sub-solar com- 
positions it is not clear that a correction is necessary. 
The above estimates suggest a composition of IM Vir 
somewhat below solar, perhaps [Fe/H] ~ —0.3, but the 
uncertainties are large and this conclusion requires con- 
firmation. 

The absolute dimensions for the system are summa- 



10 



Morales et al. 



TABLE 11 
Absolute properties of IM Vir. 



Parameter Primary Secondary 

M (Mq) 0.981 ± 0.012 0.6644 ± 0.0048 

R(Rq) 1.061 ± 0.016 0.681 ± 0.013 

T eff (K) 5570 ± 100 4250 ± 130 

logg (cgs) 4.379 ± 0.014 4.594 ± 0.017 

logL/L -0.012 ± 0.034 -0.867 ± 0.056 

BC V (mag) a -0.12 ± 0.10 -0.82 ± 0.17 

My (mag) b 4.88 ± 0.15 7.71 ± 0.29 

tisini (kms -1 ) 43 ± 2 

f S yncsini (kms" 1 )' .... 41.0 ± 0.6 26.3 ± 0.5 

a (R Q ) 5.944 ± 0.020 

Distance (pc) 89.8 ± 5.8 

a The error in these bolomctric corrections from Flowc account 
for the temperature uncertainties and include an additional 0.10 mag 
added in quadrature. 

The bolometric magnitude adopted for the Sun is M§ ol = 4.732, for 
consistency with the bolometric corrections. 

c Projected rotational velocity if synchronized with the orbital motion. 

rized in Table 1111 along with derived properties includ- 
ing the luminosities and absolute visual magnitudes. To 
calculate My we have a dopted the bolometric correc- 
tions from iFlowerl (|1996f ). with uncertainties that ac- 
count for the temperature errors as well as a contribu- 
tion of 0.10 mag added in quadrature, to be conserva- 
tive. IM Vir does not have an entry in the Hipparcos 
Catalog (|Perrvman et al.lll997t ). Here we estimate the 
distance to the system to be 89.8 ± 5.8 pc, ignoring ex- 
tinction for such a close object. Th i s is 5 0% larger than 
first estimated bv lStrassmeier et al.l (|1993f ) on the basis of 
the spectral type. Separate distance estimates for each 
component agree nearly perfectly, showing the internal 
consistency of the fundamental data on which they are 
based. The distance, proper motions from Tycho-2, and 
systemic velocity from our spectroscopic orbital solution 
(Table [2]) lead to space velocities of U = +24 .4 kms -1 , 
V = -17.0 kms -1 , and W = -3.0 kms -1 in the so- 
lar frame (with U positive toward the Galactic center). 
These do not appear to associate IM Vir with any known 
moving group in the solar neighborhood. 

With our accurate radii we compute projected syn- 
chronous rotational velocities of 41.0 ± 0.6 kms -1 and 
26.3 ±0.5 kms -1 for the primary and secondary, respec- 
tively. The value for the primary is consistent with our 
direct measurement of 43 ±2 kms -1 , suggesting that syn- 
chronization with the orbital motion has been achieved 
(if the spin axis is parallel to the orbital axis). This is 
expected from the short period. Estimates of 43 kms -1 
for the prim ary and 23 kms - 1 for t he secondary were 
reported by iStrassmeier et al.l (|1993f ). but without er- 
rors. At our request, those spectra were kindly remea- 
sured by F. Fekel (priv. comm.), giving 42 ± 2 kms -1 
and 3 1 ± 4 kms -1 , based on the calibration by iFekell 
(1997). Once again the primary agrees well with the 
synchronous value, and the secondary is probably also 
consistent, considering the difficulty of the measurement. 
An independent value fo r the primary was reported by 
IStrassmeier et all (|2000t ) as 36.2 kms -1 with an uncer- 
tainty of 2-4 kms -1 , roughly in agreement with ours. 

6. DISCUSSION 

With its solar-type primary and low-mass secondary, 
IM Vir is a particularly interesting system for testing 



models of stellar evolution. The accurate masses, radii, 
and temperatures for both components offer an oppor- 
tunity to further investigate the discrepancies for low- 
mass stars mentioned in the Introduction. The leverage 
afforded by the very different masses is unique among 
systems with at least one component under 0.8 Mq. 

We begin by comparing the measured properties 
agai nst models from the widely used Yo nsei-Yale se- 
ries (|Yi et al.ll200lt iDemarque et alj 12004 ), which treat 
convection in the standard mixing-length approximation 
with a mixing length parameter oml — 1.7432 (in units 
of the pressure scale height) , calibrated against the Sun. 
Evolutionary tracks for solar composition and for the ex- 
act masses we measure for each star are shown in the 
logg-Teff diagram of Figure flOl as dashed lines. These 
models are seen to be too hot compared to the estimated 
temperatures. Adjusting the metallicity to higher values 
brings the tracks closer to the observations, but it is not 
possible to match both stars at the same time, especially 
considering that the temperature difference is much bet- 
ter known than the absolute temperatures. Part of the 
problem has to do with the fact that these models are 
not intended for low-mass stars such as IM Vir B, which 
require a more sophisticated equation of state, and par- 
ticularly, non-gray boundary conditi ons between the in- 
terior and the photosphere (see, e.g., lChabrier &; Baraffd 
[19971) . If we therefore focus for the moment only on 
the solar-type primary, we find an excellent match to 
an evolutionary track for Z = 0.025 (corresponding to 
[Fe/H] = +0.15 in these models), at the rather old age 
of about 8 Gyr. The corresponding mass tracks for both 
stars are drawn with solid lines in Figure [TU1 and the 
shaded area indicates the uncertainty in their location 
that comes from our ~1% mass errors. An isochrone for 
this age is indicated with the dotted line. 

Not surprisingly, the secondary does not fit the evo- 
lutionary track at this metallicity, and appears implau- 
sibly old. This is a consequence of it being both too 
large and too cool compared to the models, which is in 
the same sense as deviations found for many other late- 
type dwarfs, and is believed to be due to the effects of 
chromospheric activity. Previous studies have indicated 
that better agreement for t his class of stars is possible 
with models such as those of iBaraffe et al.l (|1998l ). which 
have the required non-gray boundary conditions and use 
a lower mixing- length parameter of (Xml = 1-0. The 
latter tends to diminish (but does not completely elim- 
inate) the discrepancies noted with other models in R 
and T e ff. The one property of low-mass stars that ap- 
pears to be reasonabl y well reproduced by theory is the 
lumin osity (see, e.g., iDelfosse et all 12000 : iTorres et all 
I2002T ). This provides a means of testing the predictions 
from the Yonsei-Yale models for the primary, which point 
to an unexpected combination of super-solar metallic- 
ity and old age, and a composition that is somewhat 
inconsistent with the indications from §[5] In the left 
panels of Figure [IT] we compare the secondary lumi- 
nosity, radius , and t emperature against isochrones from 
IBaraffe et al.l (|1998| ) corresponding to the age of 8 Gyr 
found above, for different compositions including the 
value [Fe/H] = +0.15 that best fits the primary (ac- 
cording to the Yonsei-Yale models). 6 All three of logL, 

6 For this last isochrone we have extrapolated slightly in metal- 
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Fig. 10. — Evolutionary tracks for the measured masses of the 
components of IM Vir from the Yonsei-Yale series IIYi et al.H2001fc 
Dcmarquc et al. 2004). Dashed lines correspond to solar metal- 
licity, solid lines to the composition of Z = 0.025 that fits the 
primary best. The mass errors are indicated by the shaded areas. 
The implied age for the primary is 8 Gyr, and the corresponding 
isochrone is shown as the dotted line. The track for the secondary 
is calculated through the end of the main-sequence phase, which is 
reached at an age of ~40 Gyr. The age of the universe is marked 
with an asterisk on this track. The observed location of the sec- 
ondary would thus point to an implausibly old age for the star. 
This is a result of the models overestimating the temperature and 
underestimating the size for this star. 

R, and T e g for the secondary are seen to be too large 
compared to the model favored by the primary, which is 
represented by the solid lines. The measured luminos- 
ity of the secondary, as well as its effective temperature, 
suggest a significantly lower abundance, and the mea- 
sured radius is too large no matter what the metallicity. 
Changes in age at a fixed metallicity of [Fe/H] = +0.15 
do not improve the fit, as seen clearly in the right panels 
of Figure [TTJ 

There is little doubt that the radius of the secondary 
is too large compared to models. Furthermore, the 
age/metallicity inconsistency between the primary and 
secondary is perhaps an indication that the good match 
for IM Vir A in Figure [TU] is illusory, and that it too 
may have anomalies similar to those of the secondary, 
possibly related also to activity. The primary is, af- 
ter all, a rapid rotator (vs'mi = 43 kms" 1 ) . Standard 
models such as those of iBaraffe et al.l (|1998f ) do not ac- 
count for the effects of chromospheric activity, although 
the artificially low value of the mixing length param- 
eter (oml = 1. 0) seems to be a step in the right di- 
rection (see alsolTorres et alj|200l iChabrier et ai.ll20"07t 



iClausen et al.ll2009h . 

As a way of parameterizing the missing physical effects 
from activi ty and their impact on the structure of low- 
mass stars. iTorresI (|2007f ) explored the use of a correction 
factor (3 to the theoretical radii, and showed that good 
fits to empirical data could be achieved by simultane- 
ously correcting the theoretical temperatures by 




licity from the publicly available tables by Baraffc ct al. (1998) for 
[Fe/H] = 0.0 and -0.5. 



Fig. 11. — Secondary properties compared against IBaraffe et al.l 
{1993) models for parameters that best fit the pr ima ry star, when 
using the Yonsei-Yale models (see text and Figure lTOl l. Left panels: 
Fixed age of 8 Gyr, and different metallicities, as labeled. Right 
panels: Fixed abundance of [Fe/H] = +0.15, and different ages, as 
labeled. 

so as to preserve the bolometric luminosity. We apply the 
same procedure here. To account for a possible difference 
in the activity level, we have considered separate values 
of (3 for each star, and ex amined a wid e rang e of metal- 
licities and ages using the lBaraffe et al.l JT998) models to 
obtain the best simultaneous match to the properties of 
both components, at a single age and composition. In 
Figure we display the results of our grid search. A 
near-perfect fit is achieved for [Fc/H] = —0.28, as mea- 
sured by the \ 2 represented in the top panel. Interest- 
ingly, this value is much more consistent with the rough 
estimates from photometry and spectroscopy described 
in §[5] than with the metal-rich composition suggested by 
the Yonsei-Yale models for the primary. The radius cor- 
rection factors are (3\ = 1.037 and /?2 = 1.075, and the 
best-fit age is 2.4 Gyr (see middle and bottom panels). 
The age and metallicity depend almost entirely on the 
measured stellar luminosities. The values of L, R, and 
T e ff are shown in F i gure 1 131 with the adjusted isochrones 
from IBaraffe et all (|1998l ). Standard models without the 
corrections to the radii and temperatures are shown as 
dashed lines in the lower panels, for reference. 

The implication of the above fit is that the radius of 
the primary star is 3.7% larger than predicted by theory, 
and its temperature is just under 2% (100 K) too cool. 
For the secondary the models underpredict the radii by 
7.5%, and overpredict the temperature by about 3.5%, 
or 150 K. It appears, therefore, that both stars are af- 
fected by chromospheric activity. While this is not unex- 
pected for the secondary, the fact that the primary also 
seems to show the same anomalies supports other recent 
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Fig. 12. — Illustration of a grid search to obtain the best si- 
multaneous match to the mea sured properties of t he primary and 
secondary of IM Vir using the Baraffc et al. (1998) models. These 
models are adjusted here by including correction factors j3\ and /?2 
to the theoretical radii to account for the effects of chromospheric 
activity (see text). Top: The x 2 °f the fit as a function of metal- 
licity for the nominal effective temperatures of the stars (solid line 
with filled circles), and for temperatures perturbed by ±100 K to 
explore the role of systematics (dashed lines). Middle: Correction 
factors to the predicted radii. Bottom: Best-fit age as a function 
of metallicity. 

evidence that the phenomenon is not confined to low- 
mass stars, but reaches stars of solar mass (0.981 M 
for IM Vir A). Other examples of solar-type stars where 
this has been docum ented include VZ Cep B (1.108 M^; 
Torres fc Lacvll200ll ) , CV Boo B (0.968 M w :lTorres et alJ 
2008D, FL Lyr B (0.958 M^: IPopper et al 1 119861). and 
V1061 Cyg B (0.932 M^: lTorres et all 12006(1 . 

There is ample evidence of activity in IM Vir, which 
was in fact discovered through its X-ray emission, as de- 
scribed in §[T] Activity is manifested in our light curves 
in the form of spottedness, and in retrospect, our find- 
ing that the best fit to the photometry is achieved with 
one spot on each star as opposed to spots on any single 
component is consistent with the results of the stellar 
evolution modeling in the preceding paragraph, which 
indicates that both stars are affected. This is not sur- 
prising, given that both have convective envelopes and 
are rapid rotators. Additional manifes tations of activity 
in th i s case include a filled- in Ha line (jStrassmeier et al\ 
Il993t iLiu et all 119961 ) or the possible detec tion of Ha 
emission from t he secondary (iPopperl 119971 ) , Ca II H 
and K emission (iDall et al.ll2007f ). and flaring in X rays 
(jPandev k. Singhll2008f ). With the X-ray count rates and 
hardness ratio from ROSAT, t he en ergy conversion factor 
prescribed bv lSchmitt et al.1 (|1995l ). and our bolometric 
luminosities and derived distance, we have estimated the 
ratio of the X-ray to bolometric luminosity for the two 
components. We obtain log(Lx/£boi) = —3.70 ± 0.13 
and —2.84 ± 0.13, respectively, where we have assumed 
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Fig. 13.— Best-fit model from Baraffc ct al. (1998) for [Fc/H] 
= —0.28 and an age of 2.4 Gyr (solid lines), with the radii ad- 
justed by the correction factors /3i an d /32 for the primary and 
secondary, respectively (see Figure IT2l and text). The theoretical 
temperatures are adjusted by factors fj~ and the model lumi- 
nosity is unchanged. The dashed lines in the two bottom panels 
represent the standard models for the same age and metallicity, 
but without any corrections for activity (/3i,2 = 1.000). 

equal X-ray emission from each star since the ROSAT 
observation does not resolve the binary. These values 
are consistent with the secondary being completely satu- 
rated, and the primary being near saturation. They are 
similar to the X-ray luminosities seen in other active bi- 
nary systems, and are in fact in good agreement with the 
tre nd between the radiu s anomalies and Lx/iboi found 
by lLopez-MoralesI (|2007l ) (see, e.g., her Figure 2). 

The differences in radius and temperature between the 
models and the observations for IM Vir are of similar 
magnitude as those found for other active stars with 
convective envelopes, but do rely quite strongly on our 
adopted effective temperatures for the stars because of 
the way the (3 correction factors have been determined. 
Given the greater difficulty of determining absolute tem- 
peratures than other properties such as masses and radii, 
one may wonder to what extent systematic errors in T e g 
might affect the result, and also the derived age and 
metallicity. We focus here on systematic errors in the 
primary temperature, since the secondary value is tied 
to the primary through the light-curve solution. We ex- 
plored this question by perturbing the temperatures by 
±100 K, with corresponding adjustments to the lumi- 
nosities, and repeating the grid search for (3. The results 
are illustrated in the top panel of Figure Q21 where the 
dashed lines show that the best-fit metallicity changes 
by ±0.10 dex (marked by the short dotted lines). The 
changes in the (3 factors are ±0.010 for (3\ and ±0.003 
for /?2j and the change in the best-fit age is ±0.50 Gyr, 
or -20%. 
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7. CONCLUDING REMARKS 

Our extensive spectroscopic and photometric observa- 
tions of IM Vir have enabled us to determine highly ac- 
curate values for the absolute masses and radii of both 
stars to better than 2%, as well as accurate tempera- 
tures ( Table ITT|) . This eclipsing system now joins the 
ranks of those wit h the best determined properties (see 
iTorres et alJ l2009). The primary is of solar type, and the 
secondary is a late K dwarf. The very different masses 
provide increased discriminating power for testing mod- 
els of stellar evolution. We find evidence that both com- 
ponents, which have convective envelopes and are rapidly 
rotating, show discrepancies in the ir radii and tempera - 
tures compared to calculations bv iBaraffe et al. (1998), 
similar to those reported for other low-mass stars. The 
predicted radii are too small by 3.7% and 7.5% for the 
primary and secondary, and the theoretical temperatures 
are too high by 100 K and 150 K, respectively. These ef- 
fects are ascribed to chromospheric activity, for which 
there is abundant evidence in this system in the form of 
X-ray, Ca II H and K, and Ha emission, X-ray flaring, 
and spottedness seen directly in the light curves. The 
larger effect observed for the secondary is consistent with 
the larger fraction of that star being covered by spots (a 
factor of two difference). The fact that the near-solar 
mass primary is also affected supports recent findings for 
other binaries indicating that the impact of activity on 
the structure of stars is not limited to the M dwarfs, as 
previously thought, but most likely extends to all stars 
with convective envelopes, reaching up to masses near 
or larger than that of the Sun. Current stellar evolution 
models do not account for these effects. Some progress in 
this area has been made to incorporate magnetic activ- 
ity in the theoretical calculations (e.g.. iD'Antona et alJ 
12005 IMullan fc MacDonaldl[200lt IChabrier et al J 1 2007^ 
and initial results are very encouraging. 

A crucial ingredient for the comparison with models is 
the chemical composition. The best-fit models for IM Vir 
indicate an age of 2.4±0.5 Gyr and a metallicity of [Fe/H] 
= —0.28 ± 0.10. The uncertainties here represent the 
sensitivity to systematic errors of ±100 K in the effec- 
tive tem peratures, but exclu de unquantified systemat- 
ics in the Baraffc et al. (1998) models themselves, which 
are difficult to evaluate. While the metallicity prediction 
seems to be in good agreement with our rou gh spectro- 
scopic abundance estimate in §[2] and that of Pa ll et alJ 
(2007), as well as photometric estimates for both stars, 
a proper detailed spectroscopic analysis of both com- 
ponents is still lacking, and is essential to validate the 



comparison. Such determinations are notoriously diffi- 
cult in double-lined binaries, which is why relatively few 
have them. In IM Vir this task would be challenging 
for two reasons: the secondary is very faint compared 
to the primary, and it is of late spectral type. Abun- 
dance determinations for late-type stars are still prob- 
lematic due to shortcomings in the model atmospheres. 
However, IM Vir offers a unique opportunity because the 
secondary eclipse is total. High S/N ratio spectra taken 
with a sufficiently large telescope during the 28-minute 
totality phase would be of the primary only, and can be 
analyzed with standard techniques. The uncertainty in 
this observing window is estimated to be approximately 
±5 minutes, based on the uncertainties in the geometric 
light-curve parameters. 
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